Nonresonant inelastic light scattering in the Hubbard model 
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Inelastic light scattering from electrons is a symmetry-selective probe of the charge dynamics within 
correlated materials. Many measurements have been made on correlated insulators, and recent exact 
solutions in large dimensions explain a number of anomalous features found in experiments. Here 
we focus on the correlated metal, as described by the Hubbard model away from half filling. We 
can determine the B\ s Raman response and the inelastic X-ray scattering along the Brillouin zone 
diagonal exactly in the large dimensional limit. We find a number of interesting features in the light 
scattering response which should be able to be seen in correlated metals such as the heavy fermions. 
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I. INTRODUCTION 

Inelastic light scattering, and more specifically, elec- 
tronic Raman and inelastic X-ray scattering, are use- 
ful probes of the two-particle charge excitations within 
a correlated metal. Experiments are extremely difficult 
to carry out, and so far, the only correlated metals that 
have been extensively studied are the high-temperature 
superconductors 13 . An experimental challenge is to 
have a large enough energy window for the inelastic light 
scattering to be able to see the higher-energy charge- 
transfer excitations as well as low-energy particle-hole 
excitations. One class of promising materials that could 
yield interesting Raman scattering results is the heavy- 
fermion compounds that have sharply renormalized fermi 
temperatures, which compress the charge transfer exci- 
tations to a lower energy range (and thereby require a 
smaller energy window for the scattering experiments). 

On the theoretical side, much progress has been made 
in examining the inelastic light scattering of correlated 
materials. The first breakthrough occurred over a decade 
ago when Shastry and Shraiman (SS) suggested a sim- 
plification to the theory of nonresonant Raman scattering 
where in certain cases the nonresonant Raman response 
may be equal to the optical conductivity multiplied by 
the frequency. Their conjecture was recently proved 5 on 
a hypercubic lattice in d — > oo for B\ g polarization ori- 
entations. Subsequent work on the Hubbard model 6 and 
on the Falicov-Kimball model 7,5 has shown that Raman 
scattering is model independent in the insulating phase, 
and possesses all of the anomalous features seen in ex- 
periment on correlated insulators. 

In this contribution, we focus on the correlated metal, 
rather than the metal-insulator transition, and we study 
the behavior of the Raman response and the inelastic 
X-ray scattering for a wide variety of electron concen- 
trations and interaction strengths. A number of features 
arise in the Raman response that are quite interesting. 
Of course, the results for the nonresonant B\ s Raman re- 
sponse are similar to those that have been shown for the 



optical conductivity 11 ' 8 , but we discuss a number of novel 
features, including the inelastic scattering of X-rays here. 
In Section II we present our formalism, results appear in 
Section III, and we conclude in Section IV. 



II. FORMALISM 

Unlike an optical conductivity experiment, which 
probes just one symmetry of the two-particle charge ex- 
citations, inelastic light scattering can probe a number 
of different symmetries by placing polarizers on the inci- 
dent and reflected light. Three common symmetries that 
are examined are A\ & which has the full symmetry of the 
lattice, Big which is a d-wave symmetry, and £?2 g which 
is another d-wave symmetry (rotated by 45°). The Ra- 
man response that is related to the optical conductivity 
by the Shastry-Shraiman relation 4,5 is the Bi s response. 
Furthermore, the nonresonant Raman scattering in the 
i?2g channel vanishes 5 (for nearest neighbor only hop- 
ping), and in the Ai s channel it requires knowledge of 
the local irreducible charge vertex (which has never been 
calculated for the Hubbard model) , so we provide results 
only for the B\ s sector here. 

The Hubbard model 9 is the simplest model of elec- 
tronic correlations that possesses a fermi-liquid fixed 
point for a wide range of parameter space 10 . The fermi- 
liquid forms at low temperature, and often one finds quite 
anomalous behavior in the transport at higher tempera- 
tures, where the fermi-liquid coherence has not yet been 
established 11 . 

The Hubbard Hamiltonian 9 contains two terms: the 
electrons can hop between nearest neighbors [with hop- 
ping integral t* / (2Va) on a ci-dimensional hypercubic lat- 
tice 12 ], and they interact via a screened Coulomb inter- 
action U when they sit on the same site. All energies are 
measured in units of t*. The Hamiltonian is 
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where c ia (ci a ) is the creation (annihilation) operator for 
an electron at lattice site i with spin er and = c ia Ci a 
is the electron number operator. We adjust a chemical 
potential p to fix the average filling of the electrons (p = 
U/2 at half filling of p e = 1.0 for the electrons). 

Wc treat the model (1) in the infinite-dimensional limit 
so the electronic properties can be determined by dynam- 
ical mean field theory. The starting point is the locality 
of the single-particle self energy, i.e. the local Green func- 
tion is the Hilbert transformation of the noninteracting 
density of states 10 ' 8 [p(e) = cxp(— £ 2 )/y / 7t on a hypercu- 
bic lattice in d — > ool 
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with u> approaching the real axis from above. 

Formally, G(u>) can be viewed as the Green's function 
of an effective single impurity Anderson model. 10 ' 8 The 
propagator Gq(lo) for the corresponding effective non- 
interacting impurity model is then determined by 



GoV) = G- 1 (uj) + E(oj). 



(3) 



Next an Anderson impurity model solver must be em- 
ployed to determine the local Green's function from 
Gq(u) and U. Here, we use the numerical renormal- 
ization group (NRG) technique 13 ' 14 to calculate the self 
energy on the real axis, and then employ Eq. (2) to deter- 
mine the new Green's function. This algorithm is iterated 
until the Green's functions converge to a fixed point. 

The NRG is based on a logarithmic discretization 
of the energy axis, i.e. one introduces a parameter 
A > 1 and divides the energy axis into intervals 
±[A-(" +1 ),A-"] forn = 0,l,...,oo 13 > 14 . With some fur- 
ther manipulations one can map the original model onto 
a semi-infinite chain, which can be solved iteratively by 
starting from the impurity and successively adding chain 
sites. Since the coupling between two adjacent sites n 
and n + 1 vanishes like A~"/ 2 for large n, the low-energy 
states of the chain with n + 1 sites arc determined by 
a comparatively small number A sta tcs of states close to 
the ground state of the n-site system. In practice, one 
retains only these N sta ±es from the n-site chain to set up 
the Hilbert space for n + 1 sites and thus prevents the 
usual exponential growth of the Hilbert space as n in- 
creases. Eventually, after «nrg sites have been included 
in the calculation, adding another site will not change 
the spectrum significantly and one terminates the calcu- 
lation. 

It is obvious, that for any A > 1 the NRG consti- 
tutes an approximation to the system with a continuum 
of band states but becomes exact in the limit A — > 1. 
Performing this limit is, of course, not possible as one 
has to simultaneously increase the number of retained 
states to infinity. One can, however, study the A- and 
Agates-dependence of the NRG results and perform the 
limit A — > 1, A s tatos — > oo by extrapolating these data. 



Surprisingly one finds that the dependence of the NRG 
results on A as well as on the cut-off Abates is ex- 
tremely mild; in most cases, a choice of A = 2 and 
A s t a tos = 300 . . . 500 is sufficient. 

While the knowledge of the states is sufficient to calcu- 
late thermodynamic properties, the discretization leads 
to a Green's function consisting of a discrete set of poles 
and an appropriate coarse-graining procedure has to be 
applied. 15 ' 16 . Obviously, the NRG is most accurate for 
energies near the chemical potential, where the logarith- 
mic grid is finest. It is less accurate at energies far away 
from the chemical potential, because the energy grid is 
coarser and one uses an asymmetric (logarithmic) broad- 
ening for the delta functions. Since our final response 
functions involve complex convolution integrals contain- 
ing the interacting Green function, the results are ex- 
pected to be most accurate in the low-frequency region. 
At high frequencies one may still expect the qualitative 
features to be correct, but details like the precise distri- 
bution of spectral weight will typically be less reliable. 
The quantitative errors in both the Green's function and 
the inelastic light scattering response functions, however, 
are difficult to estimate. 

The nonresonant Raman response in the Bi g channel 
and the inelastic X-ray scattering along the Brillouin- 
zone diagonal (in the B\ g channel) have no vertex correc- 
tions 19 ' 7,5 and are equal to the bare bubbles. The formula 
for the Raman response has been presented elsewhere 5 
and can also be found from the SS relation 4 ' 11 ' 8 — the 
imaginary part of the nonresonant B\ s Raman response 
is 



Imxs lg (q = 0, v) = c J doj[f(uj) - f(u + v)] 

x / dep{e)A{e,ui)A{e,ui + y), (4) 



where f(oj) = 1/[1 + exp(w/T)] is the Fermi function, 
A(e,uj) = — Im[l{w + p — S(w) — e}7r] is the spectral 
function, and c is a constant. 

The inelastic light scattering by a photon with momen- 
tum q = (q,q,...,q) can also be calculated with the bare 
bubble, because one can show that all vertex corrections 
vanish here as well 20 . The result depends only on the 
parameter X = cos q, and is 
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Here we have F^z) = J dep(e)/(z — e), which is the 
Hilbert transform of the noninteracting density of states. 
Note that the X = 1 limit of Eq. (5) corresponds to 
the optical photon (q = 0) and it reduces to the Raman 
scattering case in Eq. (4) as it must. 



III. RESULTS 




FIG. 1. Evolution of the single-particle DOS for U = 4.24 
with filling for T = 0.029 (main panel) respectively with tem- 
perature for p e — 0.92 (inset). Note especially the strong 
variation of the quasi-particle peak in both cases. In addi- 
tion, with decreasing filling there is a strong redistribution of 
spectral weight from the upper Hubbard band to the Fermi 
energy. For even lower filling, the upper Hubbard band van- 
ishes completely. 

Before discussing the Raman spectra in detail, let us 
briefly look at the single-particle density of states (DOS). 
Although this quantity has been extensively discussed 
in the literature 8 it is instructive to review the major 
properties here because they will be needed to interpret 
the Raman spectra. 

The evolution of the DOS for U = 4.24 upon varying 
filling and temperature is collected in Fig. 1. Quite gen- 
erally, one can identify three structures: A comparatively 
sharp peak at the Fermi energy u = 0, and two broad 
features, the so-called Hubbard bands, representing inco- 
herent (local) charge excitations (the sharp peak vanishes 
at n = 1.0 because the system is in the Mott-insulating 
state). With decreasing filling (Fig. 1, main panel) these 
latter structures become less pronounced, the lower one 
being absorbed into the quasi-particle peak, and the spec- 
tral weight of the upper one redistributed to the Fermi 
level. For still lower filling one would eventually arrive at 



a single peak at the Fermi level, as in the non-interacting 
system. For a fixed filling, on the other hand, the varia- 
tion with temperature (Fig. 1 inset) does not affect the 
incoherent structures much, but rather leads to a strong 
(7) variation in the quasi-particle peak at the Fermi energy. 
Changing the value of U does not lead to qualitatively 
different features, only the specific values of the filling 
or temperature, where changes become pronounced, will 
vary. 

We now turn to the discussion of our results for Ra- 
man scattering by examining a correlated metal with 
U = 2.12. This system is metallic (in the paramagnetic 
phase) for all electron fillings. The fcrmi temperature is 
quite low though, so the quasi-particle peak only appears 
in the interacting DOS for low enough temperature (c.f. 
Fig. 1). This has implications for the Raman response 
as well. We show the Raman scattering at U = 2.12 for 
a variety of temperatures in Fig. 2. We plot results for 
four different electron concentrations (a) p e — 1.0, (b) 
p e = 0.9, (c) p e = 0.8, and (d) p e = 0.6. Inset in each 
panel is a blow up of the low-temperature region where 
the fermi peak appears. The behavior appears quali- 
tively different near half filling and far from half filling. 
Near half filling, we see a large charge-transfer peak at 
i/«1.5 present at high temperature, which has a fcrmi- 
liquid peak separate out of it at low temperature. The 
fermi peak narrows and is pushed to lower energy as T 
is lowered, as we expect for a correlated metal as the 
phase space for particle-hole scattering is reduced. The 
relative weight of the low energy (fermi-liquid) feature 
and the charge-transfer feature changes as we move away 
from half filling. At p e — 0.8, we find a broad high- 
energy feature, which may be separating into two peaks 
centered at approximately 1.5 and 2.5 at the lowest tem- 
perature. The scattering for p e = 0.6 is even more dra- 
matically modified. We see no large charge transfer peak 
at high temperature, just a broad low-energy peak, rem- 
iniscent of the fermi-liquid peak. As T is lowered, this 
broad peak narrows and evolves to low energy. As it does 
so, we see two additional features emerge in the Raman 
response, corresponding to a mid-infrared bump and a 
charge transfer bump. Since the interacting DOS typi- 
cally has a three peak structure (see Fig. 1), with a lower 
Hubbard band, an upper Hubbard band, and a quasi- 
particle peak, these three features arise from transitions 
within the quasiparticle peak, between the quasiparticlc 
peak and the lower band, and between the quasiparticle 
peak and the upper band. When the electron filling is 
less than half filled, the quasiparticle peak in the DOS 
lies well within what would be termed the lower Hubbard 
band. Finally, we note that the evolution of Raman scat- 
tering with doping is an evolution toward weaker correla- 
tions as p e is reduced, which is precisely what we expect 
for the Hubbard model, which is most strongly correlated 
at half filling. 

We compare these moderate correlation results to a 
more strongly correlated system (U = 4.24) in Fig. 3. 
Here three different electron fillings are shown: (a) p e = 
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1.0; (b) p e = 0.9, and (c) p e = 0.8. At half filling [panel 
(a)], the system is a correlated insulator, and the Ra- 
man response displays the characteristic features of an 
insulator, including the appearance of low-energy spec- 
tral weight at an onset temperature that is smaller than 
the T — gap, and an isosbestic point, where the Ra- 
man response is independent of temperature (here occur- 
ring near v w 3). As the system is doped, its behavior 
changes dramatically, because the doped system is metal- 
lic. Here [panels (b) and (c)] we see two peaks at high 
temperature. As T is lowered, the low-energy peak has a 
fermi-liquid peak split off and narrow and move to lower 
energy as T is reduced. The relative weight of the low- 
energy features (the fermi-liquid peak and the "mid-IR" 
peak) to the charge-transfer peak increases as the fill- 
ing moves away from half filling. We also can see that 
the fermi-liquid temperature increases as we move away 
from half filling because of the sharper fermi peak at the 
same temperature in panel (c) versus panel (b). This is 
to be expected, because as the relative correlations are 
reduced, the coherence temperature should increase. 

One might be tempted to perform an analysis similar 
to what was done for the optical conductivity 11,8 , where 
one examines the relative weights of the fermi peak, the 
mid-IR bump and the charge transfer peak. Such an 
analysis is really only feasible for the situation where 
one has a sum rule satisfied (like the optical conductiv- 
ity). The multiplication of the optical conductivity by 
the frequency to produce the Raman response obviously 
changes the sum rule, and makes such simple identifica- 
tions more difficult 21 . But the general qualitative feature 
of more scattering going into the charge transfer peak as 
the correlations increase certainly still holds. The fermi 
peak is smaller in Bi g Raman scattering, since it is con- 
centrated at small frequencies, and it is reduced by the 
extra power of frequency that multiplies the optical con- 
ductivity to produce the Raman result. 

We examine the Raman scattering at constant tem- 
perature in a moderately correlated metal U — 2.12 in 
Fig. 4. At high temperature [panel (a)], we see a broad 
charge-transfer-like feature near half filling, but as the 
electron concentration moves away from 1.0, we see spec- 
tral weight shift to low energy and evolve into a sharp 
fermi liquid peak. For intermediate fillings, we can still 
see a remnant of the charge-transfer peak at high energy, 
as a wiggle in the Raman response, but the size of that 
feature decreases sharply with filling. This picture clearly 
shows how the fermi temperature depends strongly on 
electron concentration, as the fermi temperature evolves 
to become larger than 0.341 around p e w 0.6. At low tem- 
perature [panel (b)], we see a strong fermi-liquid peak, 
that narrows and moves to lower energy with electron fill- 
ing, as expected. In addition, we see the higher-energy 
processes are sharpened, with a "mid-IR" feature occur- 
ring near v w 1 and a charge-transfer feature at higher 
energies. Once again, the relative weight of the higher- 
energy feature gets smaller as we move the filling farther 
from half filling. 



In Fig. 5, we plot the Raman response at fixed T for 
a strongly correlated system U = 4.24. In the high- 
temperature case [panel (a)], we find two prominent 
peaks, one a charge transfer peak at high energy and one 
a "mid-IR" peak at lower energy. The temperature is too 
high to see the fermi-liquid peak forming. We can see, 
however, the dramatic transfer of relative spectral weight 
to lower energy as the filling moves farther from half fill- 
ing, and the charge-transfer peak shrinks. At the lower 
temperature [panel (b)], we see a sharp charge-transfer 
peak, and an approximate isosbestic point, where the 
Raman response is almost independent of electron filling 
near v w U/2. Close to the insulating phase, there is very 
small low-energy spectral weight. As we move away from 
half filling, the low-energy spectral weight rises, devel- 
oping into a "mid-IR" feature and having an even lower- 
energy fermi peak. This result is similar to what was seen 
in the optical conductivity 11,8 , where there appeared to 
be an isosbestic point at about U/2 as a function of dop- 
ing for fillings close to half filling. The isosbestic point 
is approximately preserved here, because multiplying by 
the frequency does not remove such behavior. 

We end our results by showing the inelastic light scat- 
tering response of X-rays with momenta that run along 
the Brillouin zone boundary. Fig. 6 shows the moderately 
correlated metal and Fig. 7 shows the more strongly cor- 
related metal. Let's first start with the moderately cor- 
related results in panel 6 (a). There one sees near the 
Brillouin zone center that the fermi-liquid peak evolves 
away from the charge-transfer peak as T is lowered. But 
once the photon has finite momentum, the fermi-liquid 
feature is strongly suppressed due to the availability of 
phase-space to create particle-hole pairs. We see no peak 
separate out, but the scattering response does appear to 
sharpen as T is reduced (although it is gencrically much 
broader than what is seen at the zone center). In panel 
(b), which is even more weakly correlated, we find similar 
behavior except at the zone boundary (X = —1), where 
the results are counter-intuitive. At higher T, we see 
the low-energy X-ray scattering increase as T decreases, 
although there is only weak T dependence. But at the 
lowest temperature, the response is depleted, as if the 
system is becoming more insulating. This behavior does 
not seem to occur at other momentum transfers, or in 
the more strongly correlated case, and we do not have 
a good understanding for why it happens here, but the 
scattering of large momentum photons with correlated 
charge excitations is a complex process that is not easy 
to infer quantitative results about. 

Moving on to the more strongly correlated system in 
Fig. 7, we see some interesting behavior there as well. 
Now we can see the remnant of the fermi-liquid peak 
split-off at low energy surviving out to A = 0.5, although 
the peak does not move to lower energy as T is reduced, 
as it does at the zone center. We see the overall broad- 
ening of peaks at finite momentum, and a sharpening of 
features as T is reduced. The slope of the scattering re- 
sponse at the zone boundary can become quite large at 
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low T. 
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IV. CONCLUSIONS 

We have investigated inelastic scattering of light in 
correlated metals by solving for the dynamical response 
functions exactly in the infinite-dimensional limit. We 
examine the Hubbard model and thereby are limited to 
analyzing nonresonant scattering in the B lg sector and 
along the Brillouin-zone diagonal. The Raman scattering 
case corresponds to sitting at the zone center, while the 
inelastic X-ray scattering allows one to probe momentum 
transfers along the zone diagonal. Since the B lg response 
is related to the optical conductivity, many of our results 
are similar to those found when investigating the opti- 
cal conductivity, such as the appearance of three generic 
features corresponding to a fermi-liquid peak, a "mid- 
IR" feature, and a charge-transfer peak. The relative 
weights and temperature dependences of these three fea- 
tures vary with coupling strength and electron filling, but 
at low enough temperature, we will always see a fermi- 
liquid peak split off and evolve towards zero frequency as 
T is lowered (except at half filling with sufficiently large 
U that creates an insulator). When we allow the pho- 
ton to exchange momentum with the charge excitations, 
many of the sharp features in the Raman scattering are 
blurred, and the peaks are broadened. We no longer see 
the fermi peak evolve at finite momentum transfer, but 
as the correlations are made stronger, we do sometimes 
see remnants of the fermi-liquid peak evolution (when q 
is close to the zone center). Near the zone boundary, 
however, there is just a broad, relatively featureless re- 
sponse, whose low-energy slope increases as the system 
is more correlated and as the temperature is lowered. 

To date, there has been limited examinations of elec- 
tronic Raman scattering or inelastic X-ray scattering of 
correlated metallic systems. The most studied case is 
that of the high-temperature superconductors, but they 
cannot be examined over the full temperature range be- 
cause they have intervening magnetically ordered, or su- 
perconducting phases, which can dramatically change the 
Raman scattering. Similarly, most resonant inelastic X- 
ray scattering experiments have focused on correlated in- 
sulators, because of the difficulty in removing the elastic 
peak from the response function at the lowest energies, 
where one might expect interesting fermi-liquid features 
to appear. Our theoretical results make a number of pre- 
dictions for how the temperature, doping, and correlation 
strength dependence of inelastic light scattering varies in 
correlated metals and we believe it might be possible to 
see some of these features in heavy-fermion materials, 
which have not yet been exhaustively studied with these 
techniques. 
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FIG. 2. Raman response for different temperatures at 
U — 2.12. The different panels refer to different electron 
fillings: (a) p e = 1.0; (b) p e = 0.9; (c) p e = 0.8; and (d) 
p e = 0.6. The numbers in the legends label the temperature. 
Inset into each panel is a blow-up of the low-frequency region 
to show the fermi-liquid peak formation and evolution with 
T. 
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FIG. 3. Raman response for different temperatures at 
U = 4.24. The different panels refer to different electron fill- 
ings: (a) p e — 1.0; (b) p e — 0.9; and (c) p e = 0.8. The num- 
bers in the legends label the temperature. Inset into each 
panel is a blow-up of the low-frequency region to show the 
fermi-liquid peak formation and evolution with T. 
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FIG. 4. Raman response for different fillings at U — 2.12 
and constant temperature. The different panels refer to dif- 
ferent T, (a) T = 0.341 and (b) T = 0.0472. The numbers in 
the legends label the electron filling. Inset into each panel is a 
blow-up of the low-frequency region to show the fermi-liquid 
peak formation and evolution with electron filling. 
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FIG. 5. Raman response for different fillings at U = 4.24 
and constant temperature. The different panels refer to dif- 
ferent T, (a) T = 0.341 and (b) T = 0.0248. The numbers in 
the legends label the electron filling. Inset into each panel is a 
blow-up of the low-frequency region to show the fermi-liquid 
peak formation and evolution with electron filling. 
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FIG. 6. Inelastic 
p e — 0.8 and (b) p, 



ight scattering at U = 2.12 and (a) 
= 0.6 for different temperatures. Five 
values of the transferred photon momentum are plotted, each 
shifted by an appropriate amount, and running from the zone 
center to the zone boundary along the zone diagonal. The re- 
sults are for the Bi s sector only. In panel (a) the temperature 
decreases as the thickness decreases and ranges from 0.146 to 
0.089 to 0.054 to 0.033. In panel (b) the temperature also 
decreases as the thickness decreases and ranges from 0.240 to 
0.146 to 0.089 to 0.054. 
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FIG. 7. Inelastic light scattering at U = 4.24 and p e = 0.8 
for different temperatures. Five values of the transferred pho- 
ton momentum are plotted, each shifted by an appropriate 
amount, and running from the zone center to the zone bound- 
ary along the zone diagonal. The results are for the B\ g sector 
only. The temperature decreases with decreasing thickness 
and ranges from 0.127 to 0.077 to 0.047 to 0.029. 
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